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In this paper the implementation of data modems by means by commercially available general-purpose 


microprocessors is discussed. A number of signal-processing methods are given, which reduce and 


simplify the arithmetical operations in both the transmitter and receiver. It is shown that with these 


methods sufficient computation time can be gained to enable the use of general-purpose microprocessors 


for the design of versatile data modems for bitspeeds up to 9600 bite/s. and even higher. 


INLEIDING 


De trend naar miniaturisatie van datamodems heeft het 
onderzoek gestimuleerd naar digitale signaalbewerkings- 
technieken, geschikt voor implementatie in "custom 
made" LSI circuits of microprocessoren, voor de 
funkties van filteren, moduleren en demoduleren in 
datamodems. In het gebied van de datatransmissie 
hebben microprocessor (aP) systemen enkele specifieke 
voordelen. 

Ten eerste wordt de realisatie van "cuStom made" LSI 
circuits bemoeilijkt door de diversiteit in modulatie= 
methoden, welke door de CCITT gestandaardiseerd zijn 
voor de verschillende bitsnelheden. Deze diversiteit 
is echter geen probleem in een software georiënteerd 
AP systeem. 

Ten tweede kan een MAP oplossing economisch attrac- 
tiever zijn dan een "custom made" LSI. Het "general 
purpose" karakter van uP's maakt n‚l. massaproduktie 
mogelijk met als gevolg een snelle daling in prijs. 

Ten derde kunnen m.b.v. MAP's zeer flexibele systemen 
worden ontworpen, omdat te allen tijde programma= 
veranderingen kunnen worden aangebracht. Het grote 
aantal rekenkundige bewerkingen nodig voor de digitale 


filtering, modulatie en demodulatie overtreft echter 


al snel de mogelijkheden van de "general purpose" mP's. 


Dit maakt het noodzakelijk speciale, aan de MP aan= 
gepaste methoden te ontwikkelen voor het uitvoeren 
van de verschillende bewerkingen (speciaal de ver- 
menigvuldigingen in de digitale filters), 

Aangetoond zal worden dat deze methoden toepasbaar 
zijn voor alle gestandaardiseerde systemen voor syn- 


ehrone datatransmissie. 


SYSTEMEN VOOR SYNCHRONE DATATRANSMISSIE 


Voor synchrone datatransmissie over telefoonkanalen zijn 
een aantal modulatiemethoden gestandaardiseerd door de 
CCITT, o.a. b-fase modulatie type A of B voor 2400 bps, 
B-fase modulatie (komt in sommige opzichten overeen met 
hb # (B) modulatie) voor b800 bps en "modified" b-fase, 
keamplitude modulatie voor 9600 bps. 

Het uitgangssignaal s(t) van deze systemen met symbool 
frequentie / / / = Gy JA wordt gegeven door : 


s(é) Zu A/t- 47) sa) é „Zuhtekijengd 


wer 
waarin G). de draaggolffrequentie is, terwijl de funktie 


h(t) de spektrale omhullende van het uitgezonden signaal 
bepaalt. De relatie tussen de binaire reeks aan de 

ingang van de zender en de reeks symbolen { (xs „} 

wordt gegeven door de codeerregel welke uniek gedefinieerd 


is voor elk systeem. 


Fig. 1. Vector representatie van het gemoduleerde 


signaal op de bemonstertijdstippen. 
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De set van mogelijke waarden (x,y) van de reeks 

{Gs 7} bepaalt de bemonsterde vectorvoorstelling 
(“signal constellation") van het systeem hetgeen gra- 
fisch tweedimensionaal kan worden weergegeven zoals in 
fig. | wordt getoond voor 3 gebruikelijke systemen 


voor synchrone datatransmissie, 


METHODEN VÒÒR DE VEREENVOUDIGING VAN REKENKUNDIGE 
BEWERKINGEN IN DATASYSTEMEN 


Beschouwing van vergelijking (1) leert dat voor het 
opwekken van het lijnsignaal s(t) filtering en modu- 
latie nodig is. In een digitale implementatie ver- 
eisen deze bewerkingen (speciaal de digitale filtering 
m.b.v. niet recursieve digitale filters) een groot 
aantal vermenigvuldigingen, hetgeen een tijdrovende 
bewerking is in een AP systeem. 
Het is daarom noodzakelijk het aantal vermenigvul- 
digingen te beperken en zoveel mogelijk te vereen= 
voudigen. Een aantal alternatieven hiervoor zijn 
hieronder gegeven en zullen in meer detail worden 
toegelicht : 
1. Ontwerp van de FIR filters als interpolerende 
digitale filters. 
2. Toepassing van het principe van echomodulatie. 
3. Gebruik van "lookup tables" teneinde vermenig- 
vuldigingen in digitale filters te vermijden. 
hb, Gebruik van speciale algorithmen ("weighting 
accumulation!) voor de filters. 
5e Combinatie van een AP en een hardware vermenig- 
vuldiger. 
Sommige van deze methoden zijn algemeen toepasbaar, 
b.v. 1 en 5, terwijl de overige een meer beperkt toe- 
rassingsgebied hebben. Belangrijk is nu het verschil 
te benadrukken tussen de digitale filters in enerzijds 
de zender en anderzijds de ontvanger. In de zender 
kunnen de ingangssignalen van de filters voorgesteld 
worden d.m.v. een beperkt aantal bits, hetgeen ge- 
bruikt kan worden om de bewerkingen te vereenvoudigen, 
terwijl het ingangssignaal van de digitale filters in 
de ontvanger verkregen wordt d.m.v. A/D conversie van 
het analoge lijnsignaal s(t), hetgeen gedaan moet 
worden met tenminste 8 bits precisie. Dit verschil 
in te verwerken digitale signalen maakt het noodzake- 
lijk in zender en ontvanger verschillende technieken 
te gebruiken. Als algemeen toepasbare methode in het 
ontwerp van digitale filters kan genoemd worden het 
principe van interpolatie, waarbij gebruik gemaakt 
wordt van het verschil in bemonsterfrequentie tussen 
de ingang en uitgang van het filter om het aantal 
vermenigvuldigingen te verminderen. Indien de ver= 
houding tussen uitgangs- en ingangsbemonsterfrequentie 
b.v. gelijk is aan M en het digitale FIR filter N 


koëfficiënten heeft, dan zijn in een interpolerend 


88 


FIR filter slechts N/M vermenigvuldigingen nodig om een 
uitgangsmonster van het filter te berekenen, Is deze 
verhouding gelijk aan 1/M dan zijn M bemonsteringsinter- 
vallen beschikbaar voor het verrichten van de N vermenig= 
vuldigingen nodig voor én uitgangsmonster. Altijd moet 
getracht worden interpolerende filters toe te passen in 
zowel zender als ontvanger. 

Als tweede methode ter vermindering van het aantal ver= 
menigvuldigingen kan genoemd worden het principe van 
echomodulatie KEER waarmee de vermenigvuldigingen 
nodig voor de modulatie kunnen worden geëlimineerd. Het 
gemoduleerde signaal s(t) kan rechtstreeks worden ge- 
synthetiseerd door opeenvolgende signaalelementen te 
genereren m.b.v. een aantal parallel filters met ver- 
schillende impulseresponsies. Op deze wijze is het mo- 
gelijk modulatie te vermijden ten koste van een toename 
van het aantal bewerkingen in de filters. Dit kan voor= 
deel opleveren, als de signalen in de filters slechts 
Eén of twee bits woorden zijn. Een beperking van het 
principe van echomodulatie is echter dat het aantal 
benodigde filters slechts beperkt gehouden kan worden 
indien er een eenvoudige relatie bestaat tussen draag- 
golffrequentie @, en symboolfreguentie yy . Verder 
varieert het aantal filters en daarmee de geheugen 
capaciteit van het systeem met de verhouding a) Joop p 


De derde methode om de berekeningen in het systeem te 
vereenvoudigen is het gebruik van filters werkend met 
het principe van "look-up tables", De basis struktuur 


van dit type filter is gegeven in fig. Za. 


delay line 


a(n) | 
ebk be 


pn) 


accumulator 


p(n) 
z b 


Fig. 2. FIR filter strukturen volgens het "look-up 


table" principe. 


De inhoud van de vertragingslijn van het "finite impulse 
response! (FIR) filter wordt gebruikt als een adres 

voor een ROM, waarin alle mogelijke waarden van de 
uitgangsmonsters van het filter zijn opgeslagen. 


Op deze manier kunnen vermenigvuldigingen in het filter 


volledig vermeden worden. Teneinde de vereiste hoeveel- tellingen nodig om (2) te implementeren gereduceerd 


heid opslagcapaciteit te verminderen, welke exponentieel kunnen worden tot N simpele "tekenvermenigvuldigingen” 
toeneemt met de lengte van het adres, kunnen we het in optellingen gevolgd door slechts & volwaardige ver- 
filter in een aantal kleinere sekties opbreken (fig. 2b). menigvuldigingen met resp. A, en A, en een uiteindelijke 
De uitgangsmonsters van deze subsecties dienen nog bij optelling. Een grafische voorstelling van dit principe 
elkaar te worden opgeteld. Dit principe van “table van "weighting accumulation" wordt gegeven in fig. 3. 


look-up" worät minder attractief indien de ingangs- 
signaalmonsters in meer dan één of twee bits gecodeerd 
zijn om redenen van opslagcapaciteit en aantal vereiste 
optellingen. 

Een voobeelà van de vierde methode om de implementatie 
van FIR filters te vereenvoudigen is een speciaal ver- 
menigvuldigingsalgorithme dat we sorteerprincipe 
(engels : "weighting accumulation") genoemd hebben. 

Dit principe is gebaseerd op het feit, dat voor ieder 


van de genoemde modulatiesystemen, de ingangsmonsters 


x en y slechte een beperkt aantal aantal verschillende 
absolute waarden kunnen aannemen (bv. 0,1,3 en 5 voor 


het mod. b d / hA systeem; vergelijk fig. 1). 


Fig. 3. “Weighting Accumulator" 


Een eenvoudig pendant kent iedereen in het dagelijks 


leven. Indien iemand wil weten hoeveel geld hij in 


zijn portemonnaie heeft, terwijl de inhoud bestaat uit In een hardware realisatie lijkt dit onaantrekkelijk, 


b.v. guïdens, kwartjes en centen, dan lost men gewoon- omdat we 3 accumulatoren en additionele logica nodig 


lijk het probleem op door afzonderlijk het aantal hebben. Het principe 1s echter zeer geschikt voor AP 


guldens, kwartjes en centen te tellen en aan a Hees implementatie, waar we te maken hebben met een groot 


menigvuldigingen (met resp. 100, 25 en 1) uit te aantal registers op de chip. 


voeren plus een laatste optelling. Geven we de tijd nodig voor een vermenigvuldiging en 


In het digitale FIR filter aan de zendzijde heeft men een optelling aan met resp. T, en T, den vereist een 


te maken wet een: vergelijvsere aitueties Het uit rechtstreekse implementatie van (2) een rekentijd : 


gangssignaal P, van het filter met N coëfficiënten 


wordt gegeven door : ij Pes N ( Aa 5) 


M=/ 
Pp zi 7 zen Á, (2) Bij gebruik van het nieuwe WA algorithme en aannemende 
/ t=0 kek dat we te maken hebben met K verschillende absolute 


: : aar 1 L l 

waarin a, de ingangsmonsters voorstelt, welke TELE waarden voor de datamonsters, welke van nul verschillen 
ì tja 2 

de absolute waarden 0,1, A, en Â kunnen hebben, ter= dan wordt de vereiste rekenti) 


wijl bh, de filtercoëfficiënten met willekeurige 


waarden representeert, Ar = (Nt K) 13 % KT, (6) 


Passen we nu toe : 
omdat sgn-vermenigvuldigingen geen extra tijd kosten. 


„Nie: 
8, | : sgr ( ad Het zal duidelijk zijn dat in die gevallen waarvoor 
met : ‘-1 voor a, € 0 geldt dat K {N (zoals in de zender waar K = 3 en 
= 12 ienlij EN. liens é 
sgn{e,) = 0 voora, = Ö (3) N 0) een aanzienlijke besparing in rekentijd word 


+1 voor e, > 0 verkregen, omdat T& Et 

Betreffende de vijfde methode welke aangegeven is, 

4 N: 

wa Wiiten ve (2) herschrijven Ín de vorm, . unnen we stellen, dat het voor de hand ligt, dat een 
kombinatie van een HE en een aangepaste hardware ver= 


jh, KT n (a ) 
Pa = Ee Á, sgr (Ge zl A A, Z bi) =d menigvuldiger vele problemen kan oplossen, waarin ver= 
| a eli À menigvuldigingen het knelpunt vormen. 
n | Ts | 
De vermenigvuldiger kan dat deel van de berekeningen 
r N-# J ne /a À overnemen, welke niet goed passen in de MP met als 
t À Ee fi ig n-t gevolg een besparing van rekentijd. Deze techniek is al 
län i le Az een () toegepast in een commercieel beschikbare modem [3]- 
Hieruit blijkt dat de N vermenigvuldigingen en op= 


lake 
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ng END de PAN DALALENDENS 


De "general purpose!" AP's kunnen ruwweg verdeeld worden 
in @ groepen nl. 
— MOS AP's op een enkele chip (instructietijd = jede) 
— Bipolaire aP's \bit — sliced; instructietijd 

= 0,2 usec). 
Volgens onze ervaringen kunnen de MOS HP's (b.v. 
Intel 8080, Signetics 2650) toegepast worden voor 
datasystemen met bitsnelheden € 6800 bps, zoals h d 
en Ö f modulatie met hun eenvoudige verhouding van 
draaggolf en symboolfrequentie. 
Verder bestaan de te filteren datamonsters in deze 
systemen uit twee bits woorden, 
Hierdoor is het mogelijk echomodulatie toe te passen, 
terwijl de filters ontworpen kunnen worden als inter- 
polerende "lookup table" filters. 
In fig. b is een blokdiagram van kan JAR hd /8 gd 
systeem gegeven. 


data storage of bfilters 


filter section 
number 


interpolation 
pe number 


Fig. b. MP datazender volgens echomodulatie 


en “look-up table" filter principes. 


Het systeem bestaat uit bh "look-up table" filters 
parallel, waarvan elk weer verdeeld is in 3 kleinere 
filtersekties, De resulterende 12 elementaire filter- 
sekties maken sequentieel gebruik van één "look-up 
table" d.m.v. een multiplexer. 

Aan dit ROM worden verder 2 getallen toegevoerd, welke 
door de UP worden gegenereerd en wel een getal dat 
het filtersektie- en een getal dat het interpolatie- 
nummer aangeeft. De ME verzorgt verder de funkties 
van multiplexen en accumuleren in dit systeem. Het 
systeem bevat ongeveer 1,5 Kbytes statisch geheugen 
voor programma en tabellen, 

De uitgangsbemonsterfreguentie is 9,6 kHz. Een meer 
veelzijdig systeem kan ontworpen worden op basis van 
de bipolaire MP's (b.v. Signeties/Intel 3000 serie), 
met hun grotere arithmetrische capaciteit. Een recht- 


streekse implementatie van vergelijking (1) is gekozen 
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(fig. Sa), omdat dit een maximale systeemflexibiliteit 
garandeert. Het blokdiagram is gegeven in fig. Sb. Het 
systeem bestaat uit een differentiele coder voor de 
datasignalen, 2 identieke interpolerende digitale laag- 
doorlaatfilters, welke gebruik maken van het principe 
van “weighting accumulation" en 2 modulatoren werkend op 
een cosinus- en een sinusvormige draaggolf. Alle 
funkties van dit blokschema worden door de MHP verricht 
m.b.v. software algorithmen binnen één periode van de 
bemonsterfrequentie van 1b,b kHz. Beide filters zijn uit=- 
gevoerd met 120 filtercoëfficiënten om te kunnen voldoen 
aan de filtereisen van een modified Uh d / hb A modulatie= 
systeem (interpolatiefaktor M= 6), 

De He voert maximaal 330 micro-instructies uit om één 


uitgangsmonster van het systeem te berekenen. 


cosine 


las 
le 


Fig. 5. a} Basisschema van flexibele datazender. 
b) Blokdiagram van MP datazender met 


“weighting accumulation", 


MICROPROCESSOR IMPLEMENTATIE VAN DATAONTVANGERS 
me En ET A AN LV ANGEHG 


In de demodulator hebben we wederom te maken met de 
funkties van filtering en (de )\modulatie, 

Het is nodig de beschikbare rekentijd van het systeem te 
vergroten door toepassing van het principe van inter- 


polatie. Dit leidt tot het blokdiagram van fig. 6. 


Line 
signal data 


out 


dig. band- 
pass filter 


sine 


Fig. 6, Basisschema van flexibele demodulator, 


Post-demodulatie filters voor de onderdrukking van 
dubbele frequentie termen, veroorzaakt door het de- 
modulatieproces, zijn vermeden door gebruik te maken 
van de methode van fasecompensatie voor elk der kanalen. 
Dit maakt het mogelijk de bemonsterfrequentie in de- 
modulatoren en uitgang van de banddoorlatende filters 
te verlagen tot de Baud snelheid (interpolatie), waar= 
door de beschikbare rekentijd aanzienlijk wordt ver- 
groot. Twee digitale FIR filters met onderling go° 
fasedraaiing zijn noodzakelijk aan de ingang van het 
systeem. Deze filters moeten B-bits signaalmonsters 
verwerken, die door A/D conversie van het lijnsignaal 
worden verkregen. 

Aangezien het aantal verschillende ingangswaarden 

de filterlengte te boven gaat (M >N), kan het principe 
van "weighting accumulation" niet worden toegepast. 

Het aantal berekeningen in de filters kan echter worden 
verminderd door de twee banddoorlatende filters te 
ontwerpen door middel van simpele transformaties van 
Eén gemeenschaprpelijk laagdoorlatend filter. Noemen 

we de impulsresponsie van dit laagdoorlatend filter 

h‚ (nT) en de impulsresponsies van de bandfilters 

h‚ (nT) en h_(nT) dan kunnen we deze responsies 


schrijven als 


A, RA 
A, (15) 


A05) ola) 
Alri)sn(na,) 


) 


waarin ES de bemonsterfreguentie aan de ingang van 

de filters is (1b,h kHz) en C), de centrale frequentie 
van de filters. Kiezen we de faktor wr, gelijk 

aan 77/lh dan gaat vergelijking (7) over in : 


Á (nT‚)= A07) eos n/4|. Sqn (hr /h) 
(8) 


dh, (o75)= A loT). \ som nn/4l. Sqn (nh) 


Hieruit kan eenvoudig afgeleid worden dat, afgezien 
van de sgn=funktie, of h‚ (nT) gelijk is aan h_(nT_) 
of een van beide gelijk is aan nul. 

Derhalve kunnen de berekeningen van de van nul ver- 
schillende waarden van de impulsresponsies voor beide 
filters gemeenschappelijk verricht worden, als de 
eenvoudige sgn-vermenigvuldigingen met +1, 0 en =1 
separaat worden gedaan. De aanzienlijke reduktie van 
het aantal berekeningen, die hierdoor wordt verkregen, 
maakt het mogelijk deze flexibele demodulator m.b.v. 
een MP te implementeren. Het blokschema van het 
systeem is gegeven in fig. JT. Voor de berekening 

van een paar uitgangsmonsters (x,y) dient de MP 

2367 miero=constructies uit te voeren in 1/2b00 seconde 


in het geval van mod. kb # / hb A modulatie. 


Fig. TJ. Blokdiagram van flexibele demodulator met 
interpolerende filters en fasecompensatie 


methode. 


SLOTOPMERKINGEN 


Flexibele zender en demodulator strukturen kunnen 
worden gerealiseerd m.b.v. "general purpose! AAP!'s, 
indien aangepaste methoden voor digitale filtering en 
modulatie worden geselekteerd. In de ontvanger van 
datasystemen met hoge bitsnelheden dienen nog meer 
funkties te worden verricht zoals b.v. egalisatie. 
Op dit moment is het nog niet mogelijk zowel de 
funkties van filtering en demodulatie alsmede de 
egalisatie te doen met êên enkele MP. 
Waarschijnlijk kan dit probleem wel worden opgelost 
door gebruik te maken van de MP-vermenigvuldiger 
struktuur voorgesteld in paragraaf 3. De toename 
van rekencapaciteit kan gebruikt worden om meer 


signaalbewerkingsfunkties in het systeem te verrichten. 
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EEN EENVOUDIGE DIGITALE KORRELATIE ONTVANGER VOOR TOONSIGNALEN 


Dr.Ir. T.A.C.M. Claasen, Dr.Ir. J.B.H. Peek 


Philips Natuurkundig Laboratorium, Eindhoven 


Een digitale ontvanger wordt besproken die geschikt is voor het detecteren van toonsignalen. De struk- 


tuur van deze ontvanger is afgeleid van een korrelatie ontvanger. Er wordt echter gebruik gemaakt van 


een vereenvoudigd korrelatie principe. Het geringe aantal bewerkingen dat moet worden uitgevoerd per 


tijdseenheid maakt deze ontvanger erg geschikt voor integratie (L,.S.I.). Behalve op het principe wordt 


cok op enige simulatie resultaten ingegaan. 


INLEIDING 


Toonsignalen worden vaak toegepast voor signalerings- 
doeleinden. Voorbeelden van dergelijke toepassingen 
zijn de toondruktoets keuze bij de telefonie, oproep- 
systemen, en de zogenaamde Wmultiefrequency coding” 
(MFC) voor signalering tussen telefonie centrales. 

In al deze gevallen worden een of meerdere sinusvormige 
signalen van diverse frekwenties na elkaar of gelijk= 
tijdig uitgezonden en moet de ontvanger deze signalen 
detekteren. De frekwentie set die gebruikt wordt ligt 
vast voor elke toepassing en de ontvanger kan worden 


ingesteld voor het detekteren van een van deze frekwen= 


ties. 


Er is een uitgebreide literatuur over ontvangers voor 
toonsignalen, en ook de ontvanger die hier wordt be- 
sproken is reeds eerder beschreven KRAP Hier zal 
daarom niet al te diep op details worden ingegaan. De 
geïnteresseerde lezer wordt naar de aangehaalde arti= 


kelen verwezen. 


Vanwege diverse redenen is gekozen voor een digitale 
ontvanger. Allereerst groeit het aantal toepassingen 
waarbij de signalen in digitale vorm worden aangeboden. 
In dat geval ligt het voor de hand ook de ontvanger 
digitaal uit te voeren. Echter, ook indien analoge sig- 
malen worden aangeboden kan een digitale ontvanger 
voordelen hebben ten opzichte van een analoge ontvan- 
ger. Met name vallen te noemen de mogelijkheid van 


kostprijsreduktie door integratie (LSI) en de eenvou- 


dige programmeerbaarheid inherent aan digitale circuits. 


KORBELATIE ONTVANGER 


Uitgangspunt bij het ontwerp is de bekende korrelatie 


ontvanger voor toonsignalen waarvan een digitale versie 


is weergegeven in fig. 1. Het gedigitaliseerde ingengs= 


signaal wordt hier in twee afzonderlijke takken gekor- 


costnw,T) 


SIN( MW, T) 


Fig. 1 


releerd met respectievelijk de cosinus en de sinus van 
de te detekteren frekwentie, 

Een digitale korrelator bepaalt de som: 

N-f 

2 x(n)y(n-m) (1) 


ge) « 
n=0 
hetgeen, afgezien van een faktor 1/N, een schatting is 
van de korrelatiefunktie van de twee digitale signalen 
x(n) en y(n). 
In de ontvanger in fig. 1 worden twee van deze korrela- 
ties bepaald voor m=0, met x(n) het bemonsterde in- 
gangssignaal en y(n) gelijk aan cos(nu,T) en sin(nw T) 
respektievelijk. De zo gevonden waarden worden na kwa= 
Arateren opgeteld en de wortel uit deze som wordt dan 
vergeleken met een drempel. Detektie vindt plaats indien 
deze drempel wordt overschreden. 
Zowel de vereiste A/D omzetting als ook de vermenig- 
vuldigingen en accumulaties zijn tamelijk ingewikkeld. 
(Voor het bepalen van de wortel uit de som der kwadraten 
bestaan eenvoudige approximaties die verderop kort zul- 
len worden aangeduid). 


Er zal nu een korrelatie principe worden besproken 
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waarin deze gekompliceerde bewerkingen kunnen worden 


vereenvoudigd. 
EENVOUDIGE KORRELATIE MET HULPFUNKTIES 


Willen we een benadering maken van de funktie Pr 

uit verg. (1) dan kan dit geschieden door PE te 
maken van hulpfunkties, Uit de signalen x(n) en y(n) 
worden dan eerst signalen uln) en vln) afgeleid volgens 
de betrekkingen: 


uin) 


sgn (xens + e‚(n)) 


vn) egn ( yen) + a,(n)) 4 


waarin de Funktie sgn(x) het teken van het getal x 
voorstelt. Omdat dat dit teken alleen +1 of -1 kan 


zijn, zijn de signalen u(n) en v(n) tweewaardig (1-bit) 


De korrelatiefunktie @‚(m) is daardoor eenvoudig te 


bepalen zoals is weergegeven in fig. 2. De daarin voor- 


ain) 
bs (mi 


yin m} 


komende vermenigvuldiging kan worden verricht met een 
modulo-ë opteller en de accumulatie met een op-neer 
teller. Beide operaties zijn digitaal eenvoudig uit 


Le voeren. 


Willen we nu dat et) een benadering is van e (m) 
dan moeten natuurlijk de hulpfunkties a,(n) en a_ln) 
geschikt gekozen worden. Een mogelijke keuze voor 
(a, a) zijn realisaties van onafhankelijke witte 
ruisprocessen met uniforme amplitude verdelingsfunk- 
ties, Deze keuze is niet erg aantrekkelijk wat betreft 
de realisatie. 

Een eveneens toegelaten en aanmerkelijk eenvoudigere 
keuze zijn lineaire kombinaties van Rademacher funk- 
ties, In fig. 3a zijn 6 verschillende Rademacher 
funkties weergegeven, en het zal duidelijk zijn dat 
monsters van deze funkties eenvoudig digitaal te ge- 
gereren zijn. 

In fig. 3b zijn twee lineaire kombinaties van 
Rademacher funkties aangegeven die geschikt zijn als 


hulpfunkties voor de korrelatie, De belangrijkste 
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eisen die moeten worden vervuld zijn dat deze funkties 
geen gemeenschappelijke frekwentie komponenten mogen 
hebben met x(n) of y(n) en dat elke Rademacher funktie 
slechts in een van de hulpfunkties gebruikt mag worden. 
Aangetoond is [3] dat eu te) bij gebruik van deze hulp- 
funkties een benadering vormt van ey) afgezien van 
een konstante faktor, De benadering wordt beter naar- 
mate meer Rademacher funkties gebruikt worden in elk 
der lineaire kombinaties, In de praktijk blijkt dat 
kombinaties van drie Rademacher funkties meestal vol- 


doende zijn, 


VEREENVOUDIGDE KORRELATIE ONTVANGER 


Passen we de bovenbeschreven korrelatiemethode toe op 
de ontvanger uit fig. 1 dan krijgen we de struktuur 


die is weergegeven in fig. hb, 


cos lnw,T| 


Fig. ä# 


Omdat het ingangssignaal analoog is verondersteld, 
wordt ook een analoge hulpfunktie a, genomen. De sgn 
funktie uit vergelijking (2) wordt uitgevoerd door 


een begrenzer, waarna Éên-bits digitale monsters worden 


genomen met bemonsterperiode T. Daarbij is 2T de pe- 
riode van de hoogste orde Rademacher funktie. Op deze 
wijze is een tamelijk eenvoudige A/D omzetting ver- 
kregen. 

In plaats van de tweewaardige signalen v_{n) en v‚(n) 
te genereren op de in fig. hb aangegeven wijze kunnen 
natuurlijk ook de benodigde N waarden van elk van 
deze signalen uit een geheugen worden gehaald. Deze 
waarden kunnen immers vooraf worden bepaald omdat W 
bekend is. Dit vereist dan een 2N-bits Read Only 


Memory (ROM) of ander geheugen zoals gestippeld is aan- 


gegeven in fig. hb, 

Op deze wijze is een sterk vereenvoudigde ontvanger 
verkregen als we bovendien bedenken dat ook Zin})= 
(y £ (nh) + yin) jj op eenvoudige wijze kan worden 
benaderd door 


PACH + ly, on) | /z 
2'(n) = It Zly,enl (3) 

Iyctmsl/a + RACEN 
Iyecnsl < Lyscml- 
Verder hoeft deze waarde meestal alleen aan het einde 
van een meetperiode te worden berekend. Het zal duide- 
lijk zijn dat de zo verkregen ontvanger slechts ge- 
schikt is voor toepassingen waarbij een vaste,niet al 
te lange meettijd NT vereist is. In referentie U] 
zijn modificaties besproken waardoor de ontvanger ook 
geschikt wordt voor "kontinu bedrijf", doch hier zal 


daarop niet nader worden ingegaan. 
FREKWENTIE SELEKTIVITEIT 


Wordt aan een digitale korrelatie ontvanger als in fig. 
| die is afgestemd op frekwentie ween sinusvormig 
signaal aangeboden met frekwentie W gedurende de tijd 


NT dan zal na afloop van deze periode het uitgangssig- 


naal gelijk zijn aan 
Sin Mar 


z(N)= ESE (hi) 
2Zsin SWT 


waarin AtJ= =d, De waarde van zÍN}/N is als funktie 
van AWT uitgezet in fig. 5 voor drie verschillende 
waarden van N. Duidelijk zal zijn dat de selektiviteit 


vergroot wordt bij toenemende meettijd. 


Om de invloed van de besproken vereenvoudigingen na 
te gaan is de korrelatie ontvanger uit fig. hb gesi- 
muleerd op een digitale rekenmachine met een ingestelde 


0 Dememiibakben 


relatieve frekwentie W‚T=2E 0,1065 
zijn, opnieuw voor dezelfde drie waarden van N, weer- 


gegeven in fig. 5 


AT 


Fig. 5 


Met name voor N=5lè worden zeer goede resultaten be- 


reikt. 


CONCLUSIES 


Een digitale ontvanger voor toonsignalen is beschreven 
waarvan is aangetoond dat hij kan worden opgebouwd met 
betrekkelijk eenvoudige digitale bouwstenen. Daardoor 
is deze ontvanger erg geschikt voor integratie. Op 
details van de uitvoering van de ontvanger is niet in- 
gegaan omdat deze sterk zullen afhangen van de betref- 
fende toepassing. 

Uit simulaties is gebleken dat de ontvanger de gewenste 
werking zal vertonen mits een voldoend aantal signaal 
monsters bij de korrelatie wordt gebruikt. 

Er zijn nog een aantal vrijheidsgraden in het ontwerp 
waarmee de werking in gunstige zin kan worden beïnvloed. 
Met name de interne bewerkingssnelheid (1/T), het 
aantal gebruikte Rademacher funkties en het gebruik van 
vensterfunkties Le] zijn van invloed op het gedrag 

van de ontvanger. Nadere studie is nog vereist om de 


invloed van deze faktoren te bepalen. 
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ECHO CANCELLING IN TWEEDRAADS FULL-DUPLEX DATA TRANSMISSIE 
MET BIPOLAIRE CODERING 


Ir. H.C. van den Elzen 


Philips' Natuurkundig Laboratorium 


A mathematical description is given of an echo canceller used in a 2-wire full-duplex data transmission 


system. It is shown that for a system with bipolar encoding the advantage of optimal convergence 


properties, as is possible in a binary system,can be achieved simply by feeding the echo canceller by a 


modified binary signal rather than by the bipolar signal. Results of simulations are presented 


confirming the theory developed. An all digital implementation of the echo canceller is discussed briefly. 


INLEIDING 


Voor de overdracht van data signalen via het toekomstige 
datanet heeft de CCITT een aantal snelheden waneer). 
Een van deze snelheden is 6400 bits/sec. Als transmis- 
sieweg zullen we in dit verhaal beschouwen de verbinding 
tussen abonnee en eindeentrale (< 5 km) en de verbinding 
tussen eindcentrale en knooppunteentrale (< 25 km). In 
de laatstgenoemde verbinding kunnen Pupinspoelen voor- 
komen, zodat maximaal een spectrum tot = 3500 Hz kan 
worden overgedragen. Voor de hiervoor genoemde snelheid 
is een minimale bandbreedte van 3200 Hz nodig indien de 
bitstroom binair wordt verondersteld. Dit betekent dat 
het beschreven transmissiemedium in principe geschikt is 
voor de overdracht van 600 bits/sec in binaire vorm. 
Meestal echter worden in het transmissiemedium schei- 
dingstransformatoren toegepast zodat het gewenst is, om 
spectrale componenten in de buurt van O Hz te vermijden. 
Een codering, die "gaten" in het spectrum maakt, is de 


2) 


welbekende bipolaire codering” ‚ Kiezen we een bipolaire 
codering van de 2° orde, dan wordt de amplitude-karak- 
teristiek hiervan gegeven door [F(u)l= |sinwT|. voor 

T = 1/6400 sec betekent dit due, dat "gaten" optreden 
bij O Hz en bij veelvouden van 3200 Hz. (Bij een derge- 
lijke codering hebben we wel te maken met driewaardige 
signalen. ) 

Full-duplex transmissie wordt gewoonlijk bedreven 
over vierdraads-verbindingen. Over tweedraadsverbindin- 
gen is full-duplex transmissie evenwel ook mogelijk 
indien vorkschakelingen worden gebruikt die voldoende 
goed zijn afgesloten. Daar we echter te maken hebben 
met verschillende lengten van de verbindingen voor ver- 
schillende abonnees, is het niet mogelijk om een vol- 
doend goede compromisafsluiting van de vorken te maken. 
Dit leidt dan tot een directe lek van het gezonden 


signaal naar de eigen ontvanger en tot echos als gevolg 


van impedantie-discontinuïteiten van de kabel. We kunnen 
derhalve stellen, dat er een overdrachtsfunctie bestaat 
tussen zender en eigen ontvanger. We kunnen ook zeggen 
dat deze overdrachtsfunctie hetzg. "echo-pad" beschrijft. 
Als nu een of ander netwerk met een impulsresponsie die 
gelijk is aan de echo-pad impulsresponsie, parallel 
wordt gezet aan het echo-pad en als het uitgangssignaal 
van dit netwerk wordt afgetrokken van het ontvangen 
signaal, dan zal het gewenste ontvangen signaal (van de 
andere zijde van de verbinding) geen echo signalen, af- 
komstig van het echo-pad, meer bevatten. Het netwerk 
waar het hier om gaat wordt een echo canceller Ed 
In de praktijk kan zo'n echo canceller bestaan uit een 
transversaal filter met instelbare coëfficiënten. De 
instelling van de coëfficiënten gebeurt adaptief, d.w.z. 
het filter past zich geheel zelfstandig aan (lerend 
filter) aan zijn omgeving. 

In dit verhaal wordt een mathematische beschrijving 
gegeven van de echo canceller in een full-duplex twee- 
draads-omgeving. De beschrijving wordt gevolgd door 
enkele simulatie-resultaten en door een korte beschrij- 


ving van een digitale implementatie van de echo canceller 
SYSTEEM BESCHRIJVING 


We zullen de twee zijden van een transmissiesysteem aan 
duiden met westzijde en oostzijde respectievelijk. In 
figuur 1 is het blokschema geschetst voor de westzijde 
van een full-duplex 2-draadsverbinding met vork en echo 
canceller. 

De impulsresponsie van de zender zullen we f(t) noemen 
en het totale uitgangssignaal van de zender àlt). De 
impuls responsie van de vork vanaf de zender naar de 
ontvanger noemen we h(t), de echo-pad impuls responsie. 
Volgens figuur Ì wordt het signaal aan de ingang van de 


westontvanger : 
r(t) = e(t) - ê(t) +nalt) + n(t) (1) 
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e(t) +ns(t)+n(t) 


Fig. 1. Westzijde van een &-draads full-duplex 


transmissiesysteem met echo canceller, 


waarbij 

‚ elt) het totale echo signaal voorstelt, bestaande uit 
directe lek van zender naar ontvanger en uit echos 
als gevolg van impedantiediscontinuïteiten van de 
kabel; 

‚ &(t) het door de echo canceller nagebootste echo 
signaal 18; 

. ns(t) het gewenste signaal is, afkomstig van de 
oostzender; 


. n(t) een ruissignaal is. 


Het signaal ns(t) in (1) is het gewenste signaal voor 
de ontvanger. De taak van de echo canceller is êlt) zo 
goed mogelijk gelijk maken aan e(t), zodat r(t), afge- 
zien van de ruis n(t), zo goed mogelijk gelijk wordt 
aan ns(t), Omdat verondersteld wordt, dat de echo can= 
celler adaptief is, is het nodig dat hij informatie 
binnen krijgt omtrent het verschil e(t) — ê(t). Dit 


verschil is aanwezig in rt) en dit signaal wordt dan 


ook toegevoerd aan de regelingang van de echo canceller. 


We zien dan dat het gewenste signaal zich manifesteert 
als een "ruis" signaal voor de instelling van de echo 
canceller, De wijze waarop de echo canceller wordt in- 


gesteld komt in het volgende hoofdstuk aan de orde. 


MATHEMATISCHE BESCHRIJVING VAN DE ECHO CANCELLER 


Het datasignaal aan de uitgang van de zender (fig.1) 


kan worden geschreven als 


alt) =Ì af(t-iT), (2) 


Je 


waarin la) de datareeks met wijn +1 of -1 en T de 
bitlengte is. Nemen we nu glt) als het convolutie- 


product 
glt) = f(t) « nlt), (3) 
dan wordt het echosignaal 


e(t) = ÌÈ a,glt-jt) (4) 


J= 


De echo canceller is een transversaal filter met Nel 
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vertragingselementen, elk met vertraging T, en N 
coëfficiënten Ce De impuls responsie van de echo 


canceller is dus 


8 zi 
Rlt)= PE ec &lt-ar) (5) 
n=0. 8 


Het nagebootste echosignaal kan dan geschreven worden als 


8lt) = alt) m n(t) 


Em Eu flt-jTenT) (6) 
jee J nzo FE 

We stellen verder, dat de impuls responsie van de oost- 
zender naar de westontvanger gelijk is aan x(t). Dan 


wordt het gewenste signaal aan de westontvanger 


ns(t) = Epi, (7) 


F5 
Ds 


waarin b, de datareeks van de oostzender is (be +1 
of =1). 8 
Gewoonlijk wordt voor de minimalisatie van het 
residu echosignaal e(t) — ê(t) een mean square error 
(MSE) eriterium gebruikt op de bemonstertijdstippen 


ù 
tussen 0 en T, Vanaf dit moment zullen we de volgende 


t= tt, + mT, waarbij m een geheel getal is en to ligt 


verkorte schrijfwijze gebruiken : jen = v(t gta). Op 


grond van het bovenstaande moeten we minimaliseren 
pe ë3 a 
e= Elro} (8) 


waarbij E de mathematische verwachting voorstelt. 
Aangenomen wordt nu dat de elementen van de data reeksen 
statistisch onafhankelijk zijn en ook dat de reeksen 
la} en ,} wederkerig statistisch onafhankelijk zijn. 
De ruis wordt verondersteld een wit Gaussisch proces te 
zijn met gemiddelde OQ en variantie id (en statistisch 
onafhankelijk van de data reeksen), Gebruik makend van 
deze veronderstellingen, kan uitdrukking (Ö) met de 


hulp van (1), (bh), (6) en (T) worden geschreven als 


Mer AL 
= ES 2 Re ii 
Mi lijn mj fen Pr” 
NL! 


r 2 
Tafnj 2 find + dd (5) 


Het is zinvol om over te gaan op de volgende (getrans- 


poneerde) vector notatie 


TE = 
kl it lege: …. Je, (10) 
en li Ea dl (11) 
nn 0 | „ee N=] É ü 
met v u Ì Ps en = O(1)N-1 (12) 
nT js mit mi ® 


en de matrix notatie 
tan TED 
R = : ì (13) 
ee 
met elementen 


ka OCR 
1 = OÙ 1)N-1 


Gebruik makend van (10), (11) en (13) kan uitdrukking 


In: 


5 a Ex be} 
at A je-eFm-j-nfm-j-1? (ih, 


(9) geschreven worden als 


EE = ere = aci + & (15) 
waarbij 
2 2 2 
== : : | 
On ls Ges * A (16) 


Uit (1) volgt, dat de elementen van de matrix reëel 
zijn en symmetrisch (Hermite matrix). Gersho u) heeft 
bewezen dat een dergelijke correlatiematrix positief 
definiet is, hetgeen betekent dat de kwadratische vorm 
(15) een convexe functie is met een uniek minimum. Dit 
minimum wordt eenvoudig bepaald door de afgeleide van € 


naar c gelijk aan nul te stellen : 


grade= Ezme) =O 17) 


De optimale waarden voor c worden dan 
Se Rv (18) 
De matrix inversie die nodig is om tot de zg. Wiener 
oplossing (18) te komen is in de praktijk meestal erg 
moeilijk uit te voeren. Daarom wordt een gradient- 
methode gebruikt om iteratief de optimale coëfficiënt- 
waarden te bepalen, Als de coëfficiëntvector na k 


(xk) 


iteraties aangeduid wordt met c ‘, dan kan voor de 


coëfficiëntvector na k + 1 iteraties worden geschreven 


kk (19) 


(ki) = ct Fen 
waarbij at) een positief getal voorstelt (de iteratie- 


(k) 


factor) en grade de afgeleide is van (15) naar c na 


k iteraties. 


Voor de convergentie van c naar de optimale waarden 


k} 


zijn de iteratiefactoren ä van belang. Voor elke 


Zopt °° (x) 
iteratie is er een à' ‘zodanig, dat de totale con 
vergentietijd zo kort mogelijk wordt. De optimale 
waarden van a k) worden gevonden door (19) in te vullen 
in (15) en door de afgeleide van er naar at) 


aan nul te stellen, Het resultaat wordt dan 


TT 
grad Fi 9 
' (20) 


gelijk 


Het is duidelijk uit (20) dat voor elke iteratiestap 
een groot aantal berekeningen moeten worden uitgevoerd 
om de optimale waarden van att) te bepalen. In de 
praktijk wordt dan ook gedurende het gehele iteratie- 
proces de factor a meestal constant gehouden. De con- 
vergentie moet natuurlijk gewaarborgd blijven. Dit is in 
elk geval zo, als 


Nad = el), 0 (21) 


Met (15), (17) en (19) leidt dit dan tot 


las (22) 


graare"t) (I — aflgrade 
waarbij Ì de identiteitsmatrix is. 
Als nu de eigenwaarden van de matrix R worden aangeduid 


met A, n= O(1)N-1, dan kan uit (22) worden aangetoond 


dat voor convergentie moet gelden : 


} LI en 
a < min(z—) „n= O(1)N-1 (23) 
n 
De convergentiesnelheid is uiteraard niet optimaal, daar 
geen optimale waarden voor à worden gebruikt. Voor wat 


k) 


betreft deze convergentiesnelheid heeft Gersho aange- 


toond, dat een goede praktische waarde voor Q is 


a = (A ei sj (ah) 

max mir 

met À _ en À . de maximale resp. de minimale eigen- 
max min 


waarde van de matrix È. 
DE CORRELATIE MATRIX R 


We bekijken nu eerst hoe de matrix E eruit zou zien, 

als we te maken hadden met een binair systeem. De impuls 
responsie van de zender is dan ê(t) en met de hulp van 
(Ih) vinden we dan R= I. In dit geval is het dan ook 
gemakkelijk om de optimale waarde voor a te vinden uit 
(20); Base * 5. 


Uit (17) en (18) volgt verder : grad Ee = ale-v) en 


Ek * E Toepassing van (19) levert : 


NSE, u be Balen) = V (as) 


Dus na Één iteratie worât, in dit geval, reeds de 
optimum coëfficiëntvector gevonden, onafhankelijk van de 
beginwaarde van deze vector, 

Bij de bipolaire codering van orde 2, wordt de 


impuls responsie van het lineaire gedeelte van de zender, 
f(t) = êlt) - &{t-2T) (26) 


Met (26) en (1b) vinden we dan eenvoudig de elementen van 


de matrix R : 


5 ek E En „1+2 De Stal (ar) 
waarbij 
Tei = d 
ij loi 4 E ‚ is de Kroncker delta. 


De matrix R wordt derhalve : 


2 =| 0 Dee 0 

| 

0 Gt: 8 f 

1 SD Oel ' 

LI 

Ral Ol 0 2 0 & (28) 

oe lt 0 8 

LI LI 

1 1 

2 


Q-eee 


De eigenwaarden van zo'n matrix werden afgeleid door 
Diane”. Het resultaat is, dat de som van de maximum en 
minimum eigenwaarde ongeveer h is. Volgens (2l) hebben 

we dan q= Ì. Zoals reeds gezegd is deze a niet optimaal 
en er zullen meerdere iteraties nodig zijn om de optimale 
waarde van de coëfficiëntvectorte vinden. De convergentie. 
snelheid voor een bipolair systeem is dus kleiner dan 
voor een binair systeem, als beide systemen volgens fig] 


zijn opgezet. 
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PRAKTISCHE INTERPRETATIE 


Het zendgedeelte van een systeem met bipolaire codering 


bestaat uit een binaire modulo-2 transformatie om fouten. 


propagatie in het ontvanggedeelte te voorkomen, en een 
lineaire operatie voor de spectrumbeheersing. Als we 
het convergentievoordeel van een binair systeem be- 
schouwen, kunnen we ons afvragen of het misschien mo= 
gelijk is om dit voordeel ook te behalen bij een 
bipolair systeem, Het antwoord is ja, immers, als we het 
binaire signaal na de modulo-2 transformatie aan de echo 
canceller toevoeren, i.p.v. het ternaire signaal na de 
lineaire operatie, dan kan in alle formules vanaf (5) 

de impuls responsie f(t) vervangen worden door 8(t) en 


daardoor verkrijgen we voor R weer de identiteitsmatrix. 


Fig. 2. Westzijde van een Z-draads full-duplex 
bipolair gecodeerd data transmissiesysteem 


met echo canceller. 


Dus, hoewel we bezig zijn met een bipolair transmissie- 
systeem, gebruiken we toch het convergentievoordeel van 
een binair systeem. In figuur 2 is weer de westzijde 
weergegeven van een full-duplex 2-draads systeem met 
bipolaire codering. Met de schakelaar in positie II 
wordt aan de echo canceller het bipolaire signaal toe- 
gevoerd (we zullen dit geval "Bip cancel" noemen), zo- 
dat (28) geldt. Met de schakelaar in positie TI wordt 
aan de echo canceller het getransformeerde binaire sig=- 
naal toegevoerd (dit geval zal "canceller" genoemd 
worden). Hier hebben we dus te maken met de identiteits 
matrix, met de mogelijkheid tot optimale a. 

In de mathematische beschrijving werd de gradient- 
vector exact bekend verondersteld, In de praktijk 
echter is dit niet mogelijk, aangezien de "mathema- 


tische verwachting" benaderd dient te worden door een 


100 


tijdgemiddelde, zodat (8) verandert in 


Ü K ë 
E mn 


met K een (groot) getal dat aangeeft hoeveel bits ge- 


bruikt worden voor de middeling, en 


> A (30) 


De afgeleide van € naar de i° 
wordt dan met (29) en (30) 


A eo 
graded AE Esfage 00 
é Jeen 


De laatste uitdrukking zegt ons, dat een benadering van 


coëfficiënt Cs i=0( 1ÌN-1 


de gradient wordt verkregen door het ontvangen signaal 
ee te correleren met de vertraagde versies van het 
ingangssignaal aan de echo canceller. Het blokschema van 


de echo canceller wordt dan zoals in figuur 3 aangegeven. 


Fig. 3. Basisstructuur van de echo canceller, 


VOORBEELD MET EEN COEFFICIENT 


Een eenvoudig voorbeeld met slechts êên coefficient moge 
de voorgaande theorie verduidelijken. 


We nemen (overeenkomstig (15)) : 
e= 5 - he + 5 (32) 


Het minimum van e wordt gevonden via S£ = 2e - he 0; 
opwel Sopt” 2 en Ein = |. 

Met de gradient-methode wordt het optimum van de coëf- 
ficiënt iteratief gezocht volgens (19) : 


ED U) _ 0) roolk)_ 1) (33) 


Substitutie van (33) in (32) levert 
elk) el), atnid (de tP) 4)? (31) 


n Nr + 
Uit (3l) is eenvoudig te vinden dat eik 1 elk) 


<0, als 
O< asl, met dS } (vergelijk ook (20) en (23)). Inde 
figuren hb en 5 zijn slaap. sk uitgezet tegen het 
aantal iteraties k met als parameter de iteratiefactor 

a (de punten zijn onderling verbonden door rechte stip= 


pellijnen). 
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Fig. hb. Het verloop van de coëfficiëntwaarde tegen 


het aantal iteraties met a als parameter. 
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Fig. 5. Het verloop van € tegen het aantal itera- 
ties met a als parameter. 


Om na te gaan wat er gebeurt als de gradient niet exact 


Oe) 


gelijk genomen wordt aan (2c werd een verstoring 


aangebracht na de eerste iteratie : 


ij (1) 


= grade ‘+1, In figuur 6 is aangegeven hoe 


(grade 

deze verstoring uitwerkt op het verloop van elk), We 
constateren, dat door het nemen van een verkeerde gra- 
dient een tijdelijke divergentie optreedt in het verloop 
van E « We zien ook, dat de verstoring ernstiger is naar 
mate a groter is (vergelijk bv. het verloop voorad= “en 
a= ff in fig. 5 en fig. 6). Doordat in een praktische 
realisatie de gradient verkregen wordt door (eindig) 
correleren en derhalve slechts een benadering van de 


werkelijke gradient is, zullen praktijkcurves een soort- 


gelijk karakter hebben als die in fig. 6. 


Fig. 6. Het verloop van € tegen het aantal iteraties 
met a als parameter. Na één iteratie wordt 


een fout in de gradient verondersteld. 


SIMULATIE-RESULTATEN 


Een typische impuls responsie x(t) voor een systeem met 
bipolaire codering is gegeven in figuur 7. De totale echo 
impuls responsie glt) is eveneens in de figuur weergege- 
ven. Voor de simulaties werd als bemonstertijdstip 


to 5 + mT genomen. Het aantal coëfficiënten van de echo 


Fig. T. Impuls responsies van : 
. het gewenste signaal (x(t)); 


. het ongewenste echo signaal (g(t)). 


canceller was N= 16 en de ruis werd verondersteld niet 
aanwezig te zijn. Gebruik makend van (19) en (31) werd E 
berekend volgens (29) met a als parameter. Voor de data- 
reeksen la} en (b} werden verschillende pseudo-random 
reeksen gebruikt. 

Figuur B geeft enkele (genormaliseerde) resultaten weer 
in het geval van simplex transmissie, ofwel ns(t)=0; al- 
leen de westzijde zendt. Voor de bepaling van de gra= 
dient werd gecorreleerd over 300 bitintervals. We zien, 
dat voor "Canceller" (positie I van de schakelaar in 
fig.2) met a = 5 na Één iteratie het minimum van € nage- 


noeg 1s bereikt, d.w.z, een correlatietijd over 300 bits 
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ner, In figuur 9 worden enkele (genormaliseerde) resul- 
dn taten getoond voor full-duplex transmissie met een cor= 
relatietijd over 2400 bitintervals. Voor "Canceller" met 
a = 3 daalt e na een iteratie tot een bepaald niveau, 
doch volgende iteraties geven geen verdere daling te zien. 
Dit betekent, dat een correlatietijd over 2400 bits hier 
nog niet voldoende is, Een langere correlatietijd is ge- 


wenst, ôf een kleinere a ‚zoals bv. a = Ì. Voor "Bipcan- 


1 1 
0 2 | cel" werd weer genomen a = B FeSD.5 
ee k 


a bie Uit de figuren 8 en 9 blijkt, dat de convergentie- 


Bipcancel snelheid voor "Canceller" altijd beter is dan voor "Bip- 


cancel", hetgeen een bevestiging is van de theorie. 


EEN DIGITALE IMPLEMENTATIE VAN DE ECHO CANCELLER 


In figuur 10 wordt een blokschema getoond van een digi- 
taal uitgevoerde echo canceller volgens het in dit verhaal 
beschreven principe. Omdat volgens figuur 3 een groot 


aantal correlatoren nodig is, werd gekozen voor een 


Fig. 8. Simplex transmissie : genormaliseerde mean 
square error tegen het aantal iteraties 


met 300 bitintervals per iteratie. 


is in het simplex geval voldoende lang om een goede 
schatting van de gradient te maken. Voor a = ù conver- 
geert het systeem langzamer naar zijn eindwaarde, zoals 
te verwachten was. 

Simulaties met a = à voor "Bipcancel" (positie II van de 
schakelaar in fig. 2) resulteerden bij een correlatie- 
tijd over 300 bits in een verloop dat soms nog tijdelijk 


divergentie vertoonde, alhoewel het geheel convergeerde. 


Voor “Bipcancel" moet daarom de correlatietijd langer 


worden genomen of zoals in fig. Ö de waarde van à klei 


E” er Fig. 10. Digitale implementatie van de echo can= 
{ a=% at=é 
Canceller 20 Bipcancel celler. 


sequentiële structuur, zodat slechts Één vermenigvuldiger 
en accumulator nodig is, die voor alle correlatoren wor= 
den gebruikt. Voor N coëfficiënten en derhalve N correla- 
toren betekent dit dat de interne snelheid van het appa- 
raat N keer hoger is dan nodig zou zijn voor de parallel- 
structuur volgens fig. 3. In elk bitinterval T komt een 
nieuwe data-sample binnen in het dataregister. Dit re- 
gister is uitgevoerd als een recirculerend register met 
lengte N en een schuiffregquentie N/T. De coëfficiënten 
schuiven met dezelfde frequentie in het coëfficiëntenre- 
gister, dat eveneens een recirculerend register is. Als 
de coëfficiënten worden opgeslagen als M bit seriewoorden 
is de interne schuiffrequentie van het coëfficiëntenre- 
gister gelijk aan M.N/T. In één bitinterval T verschijnen 
er dus N data-samples en N coêfficiënten aan de uitgangen 
van de registers, Ze worden met elkaar vermenigvuldigd en 
Fig. 9. Full-duplex transmissie: genormaliseerde daarna opgeteld in de accumulator. Aan het einde van elk 


mean square error tegen het aantal itera- bitinterval is er zodoende een echo-sample benaderd. De 


ties met 2400 bitintervals per iteratie. benaderde echo-samples worden afgetrokken van de ontvan 
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gen signaal-samples, die verkregen worden uit het ont- 
vangen signaal via een A/D omzetter (bemonsterfrequentie 
5. Het uitgangssignaal van de aftrekker wordt naar de 
correlator gevoerd (Één waarde per bitinterval). De 
andere ingang van de correlator is het uitgangssignaal 
van het dataregister (N samples per bitinterval). De 
correlator vermenigvuldigt beide signalen en telt de pro- 
ducten op bij de reeds aanwezige waarden, Na K bitinter- 
vals worät het uitgangssignaal van het correlatorregis- 
ter vermenigvuldigd met a en het (K+1)® bitinterval 
wordt gebruikt om de N coëfficiënten in het coëfficiën- 
tenregister te veranderen volgens (19) en om het corre- 
latorregister weer op nul te zetten. 

De echo canceller is gebouwd met N = 32 coëfficiënten. 
De interne verwerking geschiedt met M= 16 bits, wat 
leidt tot een snelheid van 3,2 MHz voor 6h00 bits/sec. 
De correlatietijd is instelbaar voor 28, 2), 21 of al 
bitintervals, terwijl de factor a instelbaar is voor de 


waarden zl, gs En 


CONCLUSIES 


In dit verhaal werd een data transmissiesysteem met bi- 
polaire codering behandeld, waarmee het mogelijk is om 
6400 bits/sec over te äragen via gepupiniseerde verbin- 
dingen. Om full-duplex te kunnen werken over € draden, 
werd het systeem voorzien van een echo canceller, waar= 
van een mathematische beschrijving werd gegeven. Met 
deze beschrijving was het mogelijk om te bewijzen dat, 
als de echo canceller gevoed wordt met het bipolaire 
signaal, de convergentiesnelheid altijd lager is dan zij 
zou zijn in een binair systeem. Als echter de echo can= 
celler gevoed wordt door het binaire signaal direct na 
de niet-lineaire transformatie die de bipolaire codering 
voorafgaat, dan kunnen dezelfde convergentieeigenschap- 
pen worden behaald als in een binair systeem. 

De theorie werd bevestigd door computersimulaties en 
middels een praktisch gerealiseerd model van de echo 


canceller in digitale vorm. 
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EEN OVERZICHT VAN NIET-LINEAIRE EFFECTEN IN REKURSIEVE DIGITALE FILTERS 


Dr.Ir. T.A.C,.M. Claasen, Dr.-Ing. W.F.G. Mecklenbräuker, Dr.ir. J.B.H,Peek 


Fhilips Natuurkundig Laboratorium, Eindhoven 


Een overzicht wordt gegeven van kwantisatie en overloop effecten in rekursieve digitale filters die 


werken met vaste komma getallen representatie. Een aantal structuren die zijn voorgesteld voor het voor- 


komen of reduceren van limiet-cycli en overloop verschijnselen worden besproken, 


INLEIDING 


Bij de realisatie van elk rekursief digitaal filter is 
het noodzakelijk twee verschillende niet-lineaire be- 
werkingen uit te voeren. 

1) Kwantisatie, die moet worden uitgevoerd om de woord- 
lengte van de signalen die is toegenomen na bijvoor- 
beeld een vermenigvuldiging weer op de oude waarde 
terug te brengen. 

2) Overloop nietlineariteiten die de signalen zodanig 
moeten beperken dat ze in absolute waarde nooit 
groter kunnen worden dan het maximale getal dat kan 
worden voorgesteld in de gekozen getallen represen- 
tatie. 

Beide niet-lineaire bewerkingen zijn noodzakelijk doch 

introduceren ongewenste effecten in het digitale filter. 

Deze effecten zijn uitvoerig beschreven in de litera= 

tuur [',2,3). en enkele ervan zullen hier nader worden 

toegelicht. Verder zal worden ingegaan op enige me- 
thoden die zijn ontwikkeld om deze ongewenste verschijn= 
selen te verminderen of te elimineren. 

In dit artikel zullen we ons voornamelijk concentreren 

op een tweede orde rekursief digitaal filter, en de 

zogenaamde "wave digital filters (WDF)". WDF's zijn 
reeds in de door Prof. Fettweis gegeven presentatie 

aan de orde geweest. Het tweede orde filter (2e orde 

sektie) is een standaard bouwblok waaruit hogere 

orde filters kunnen worden opgebouwd, en is daarom erg 


belangrijk. 


Er zijn diverse structuren bekend voor de opbouw van 
deze 2e orde sektie. Het rekursieve deel van een veel 
gebruikte structuur is dat uit fig. 1. 

In deze figuur zijn de mogelijke plaatsen voor de 
kwantisatoren (Q) en de overloop nietlineariteiten (P) 
aangegeven. Voor de nietlineaire effecten die hier 
zullen worden besproken zijn alleen de nietlineari= 


teiten van belang die in het rekursieve deel (terug- 


Tijdschrift van het Nederlands Elektronica- en Radiogenootschap deel 42 -— nr. à& — 1977 


koppel lussen) aanwezig zijn. 

Voor de plaatsing van de kwantisatoren zijn er dan 
twee verschillende mogelijkheden: a) direkt na de ver- 
menigvuldigingen doch voor de optelling (azlalex), 

b) na de optelling (CQ, (x)=Q, Cadex). 

Verder kan de kwantisatie nog op diverse manieren wor= 
den uitgevoerd. Twee van deze manieren zijn weer- 
gegeven in fig. 2. Bij afronden (fig. 2a) wordt het 
dichtsbijzijnde getal genomen, bij afhakken (fig. 2b) 
worden de overtollige bits eenvoudig weggelaten. 

Na optelling van de getallen ay(n-i), byln-2}) en uln) 
is het mogelijk dat deze som groter is geworden dan 
het maximale woord dat in het register past. In dat 
geval wordt door de nietlineariteit P een ander getal 
afgeleverd dat wel in het register past. Ook hiervoor 
zijn er weer diverse mogelijkheden, waarvan er drie 
zijn geschetst in fig. 3. 

Het eenvoudigste is natuurlijk om de overloop veroor= 
zakende bits maar weg te laten (fig. 3c) of het 
signaal bij overloop nul te stellen (fig. 3b). Het 
zal blijken dat deze oplossingen aanleiding kunnen 
geven tot niet gewenste verschijnselen. Daarom zijn 


andere oplossingen voorgesteld zoals de verzadigings 
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nietlineariteit uit fig. 3a. 
LIMIET CYCLI 


Limiet cyclill.c.) zijn periodieke trillingen die kun- 
nen ontstaan in digitale filters als gevolg van de 
niet-lineaire bewerkingen in het rekursieve deel bij 
afwezigheid van een ingangssignaal (x(n}z0). Te onder- 
scheiden zijn l.c. die worden veroorzaakt door de 
kwantisatie en l.c. afkomstig van de overloop. Deze 
laatste worden meestal overloop oscillaties genoemd 

en zullen later aan de orde komen, 

In fig. hb zijn 16 verschillende l.c. weergegeven die 
zijn gevonden in een Ze orde sectie als in fig. 1, met 
Eén afrondingskwantisator. Ze zijn gevonden door het 
filter telkens te starten met andere beginvoorwaarden 
(inhouden van de registers). De grootste l.c, heeft 


een amplitude van 90 kwantisatie stappen. 


Fig. } 
Het zal duidelijk zijn dat deze oscillaties erg storend 
kunnen zijn in een digitaal filter. Immers bij pauzes 
in bijvoorbeeld spraak signalen zal het filter spontaan 
gaan oscilleren. 
Aangetoond kan worden dat bij gebruik van afronden dit 
soort oscillaties altijd zal optreden in bijna elk 
filter van praktisch belang. Er dienen dus maatregelen 
genomen te worden om deze l.e. te onderdrukken of te 
vermijden, waarvan er nu enkele worden besproken. 
1) Voor de hand ligt om de kwantisatiestap grootte q 
te verkleinen waardoor de amplitude van de l.c. zal 
afnemen. Dit kan gebeuren door meer bits te nemen 
voor de signaal-representatie, Enerzijds verhelpt 
dit natuurlijk niet de kwaal, l.c. blijven bestaan, 
anderzijds zijn vaak een groot aantal bits nodig 
om de l.c. amplitudes zo klein te krijgen dat ze 
niet meer storend zijn. Een groot aantal bits ver- 
kleint de rekensnelheid en kompliceert de realisatie, 
zodat andere oplossingen gewenst zijn. 
2) Omdat we hebben gezien dat afronden altijd tot l.c. 
leidt is het van belang te kijken wat er gebeurt als 
de signalen worden afgehakt (fig. 2b). Worden twee 


van deze kwantisatoren gebruikt dan blijkt dat alleen 


3) 


zeer specifieke l.c. in de 2e orde sektie kunnen op- 
treden namelijk l.c. met periode 1 (konstant uit- 
gangssignaal) en periode 2 (alternerend uitgangs- 
signaal). Deze treden alleen op in laag en hoog 
doorlaat filters, terwijl filters met een doorlaat- 
band rond een kwart van de bemonsterfrekwentie l.c. 
vrij zijn. Bouwen we nu een hogere orde filter op 
door meerdere van deze sekties in cascade te plaat- 
sen dan zullen l.c. die bijvoorbeeld in de eerste 
sektie zijn opgewekt volgende sekties aanstoten 
waardoor opnieuw l.c. van diverse amplitudes en 
perioden aan de uitgang van het filter kunnen op- 
treden. 

Een geheel andere situatie treedt op als slechts 
een kwantisator wordt gebruikt waarin afhakken 
plaats vindt. Voor bijna alle voorkomende waarden 
van de filter coëfficiënten zal de 2e orde sektie 
vrij zijn van l.c. Dit is weergegeven in fig. 5 


waarin een äriehoekig gebied is aangegeven.in het 


„ seh meantinke 
ta Mms 


fd 


U PS 


Fig. 5 


vlak van de coëfficiënt-waarden a en b. Alle voor= 
komende filters hebben coëfficiënt waarden die lig- 
gen binnen deze driehoek. Limiet eycli zullen alleen 
optreden voor sommige laag- (en hoog-) doorlaat 
filters met hoge kwaliteitsfactoren. Voor deze 
filters liggen de coëfficiënten in de gearceerde 
trapeziumvormige gebieden, waarvan er een vergroot 
is weergegeven. De coëfficiënt-waarden zijn daarbij 
voorgesteld met 8 bits en aangegeven zijn het aantal 
mogelijke l.c. voor elk coëfficiëntenpaar (a,b). 
Waar niets staat aangegeven zijn ook geen l.c. ge- 
vonden. 

Een andere methode om van l.c, af te komen is de 
kwantisatie te laten besturen door een hulpsignaal 
dat ofwel uit in het filter aanwezige signalen 
wordt gevormd, ofwel aan een (pseudo-) toevals 
generator wordt ontleend. In het eerste geval moet 
het stuursignaal zorgvuldig worden gekozen zodat 
inderdaad geen l.c. meer optreden. In het tweede ge- 
val zorgt het toevalskarakter van het stuursignaal 
ervoor dat de regelmaat van de l.c. wordt verstoord 


en het resulterend gedrag ruisachtig van aard is, 


hetgeen over het algemeen minder storend is. Beide 
maatregelen vereisen extra “hardware en leiden daar- 


om tot een meer gekompliceerde realisatie. 


Tot nu toe hebben we alleen de 2e orde sektie beschouwd. 
Ook in WDF's kunnen echter l.c. optreden indien de sig- 
nalen worden afgerond. Er is aangetoond dat in deze 
filters l.c. altijd kunnen worden vermeden indien de 
overtollige bits van de signaalmonsters worden afge- 
hakt juist voor ze in de registers worden geplaatst. 
Gebruikmakend van dit feit zijn 2e orde sekties uit 
WDF's afgeleid die eveneens l.c. vrij zijn bij gebruik 
van afhak kwantisatoren. Een voorbeeld van een derge- 
lijke structuur zal door Prof. Butterweck en Ir. Ver- 
kroost nader worden toegelicht. Deze sektie heeft als 
additioneel voordeel dat l.c. zelfs kunnen worden voor- 
komen indien een konstant ingangssignaal aan het filter 
wordt aangeboden, doch dan moet de kwantisatie gestuurd 


worden door het ingangssignaal, 
OVERLOOP OSCILLATIES 


Zoals reeds vroeger is opgemerkt kunnen ook de overloop 
nietlineariteiten aanleiding geven tot oscillaties 

bij afwezigheid van een ingangssignaal. Een voorbeeld 
van een dergelijke oscillatie is weergegeven in fig. 6 
waarin dezelfde oscillatie driemaal is weergegeven, 
telkens met een andere tijdsschaal. Deze oscillatie 

is gevonden in een 2e orde sektie die was gerealiseerd 
met een overloop karakteristiek als in fig. 3b is 
aangegeven. Omdat de amplitude van dergelijke oscilla= 
ties in de buurt ligt van het overloop niveau, zijn 
deze oscillaties natuurlijk te allen tijde ontoelaat- 
baar. Ook bij gebruik van een overloop nietlineariteit 
uit fig. 3e zullen zulke oscillaties kunnen optreden, 
zodat in de praktijk deze eenvoudig te realizeren 
overloop nietlineariteiten niet kunnen worden toegepast. 
Aangetoond is dat overloop oscillaties niet zullen 
optreden in een 2e orde sektie als de karakteristiek 
van de overloop nietlineariteit in het in fig. T aan- 
gegeven gebied ligt. 

Duidelijk zal zijn dat de in fig. 3a weergegeven ver- 
zadigingskarakteristiek hieraan voldoet en bovendien 
net voordeel heeft dat van alle toelaatbare karakteris- 
tieken deze de kleinste fouten geeft. 

Voor WDF's is aangetoond dat overloop oscillaties 
nooit zullen optreden onafhankelijk van de gebruikte 


overloop karakteristiek. 
ANDERE OVERLOOP PHENOMENEN 
Tot nu toe is alleen het geval beschouwd dat geen in- 


gangssignaal aanwezig is. Met ingangssignalen wordt 


de situatie echter anders. 
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Ter illustratie van de verschijnselen die nu kunnen 
optreden gaan we uit van een 2e orde sektie met ver- 
zadigings nietlineariteit. Deze sektie wordt geëxci- 
teerd met een sinusvormig signaal met amplitude A en 
frekwentie f. 

We laten nu eerst langzaam (relatief t.o.v. de periode) 
de amplitude van het ingangssignaal toenemen. Het 
uitgangssignaal zal eveneens sinusvormig zijn, en ook 
de amplitude daarvan zal dus langzaam toenemen. Dit is 
weergegeven in fig, Öa, waarin een sterk gekompri- 
meerde tijdsschaal is gebruikt. Bij voldoend grote 


uitsturing zal het uitgangssignaal vastlopen door de 
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aanwezigheid van de verzadigingsnietlineariteit. Dit 
is duidelijk te zien in fig. Ba. Vanaf dat moment 
wordt de waarde van A weer langzaam verminderd, We zien 
dan echter dat de amplitude van het uitgangssignaal 

op de maximale waarde blijft hangen. Dit blijft zo tot 
dat de ingangsamplitude beneden een bepaalde waarde 
komt, Op dat moment springt de uitgangsamplitude plot- 
seling naar de waarde die cok werd aangenomen bij de- 
zelfde ingangsamplitude maar voor het vastlopen. Dit 
effect wordt daarom in de engelstalige literatuur met 
“jump phenomenon" aangeduid. Na de sprong volgt de 


uitgangsamplitude opnieuw de ingangsamplitude. 


ingang en 


uitgang 


| a) 


amplitude 
aan uitgang a\ 


—» frekwentie 
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Een zelfde sprongverschijnsel kunnen we ook konsta= 
teren indien bij konstante amplitude A de frekwentie 

f langzaam wordt gevarieerd. In fig. 8b zijn twee ge- 
vallen weergegeven. In geval a was A zodanig dat het 
uitgangssignaal ook bij de resonantiefrekwentie niet 
vastliep, en a geeft dus de resonantiekromme van de 

ee orde sektie weer, In geval b was A groter genomen en 
liep het uitgangssignaal wel vast, Ook hier bleef het 
uitgangssignaal "hangen" op het maximale niveau en trad 
een plotselinge sprong op bij verder toenemende fre- 
kwentie f. 

In de 2e orde sektie zullen deze effecten kunnen op= 


treden ongeacht de gebruikte overloop karakteristiek, 


Het minste last geeft de verzadigingskarakteristiek, Press, New York, 1975. 


waarvoor in fig. 9 de waarden van de coëfficiënten 2 J.F. Kaiser, "On the limit cycle problem", in Proc. 
a en b zijn weergegeven waarvoor deze effecten niet 1976 IEEE Int. Symp. Circuits and Systems, Munich, 
optreden. Germany, 1976, pp. 6h2-645, 


3 T.A.C.M., Claasen, W.F.G. Mecklenbräuker, J.B.H. 
Peek, "Effects of Quantization and Overflow in 
Recursive Digital Filters", IEEE Trans. on Acous- 
tics, Speech, and Signal Processing, Vol. ASSP- 
eh, Ho. 6, December 1976, pp. 517-529. 


Fig. 9 


Speciaal hoog- en laagdoorlaat filters hebben erg veel 
last van deze verschijnselen. Voor deze filters zijn 
deze verschijnselen te onderdrukken door toepassing 
van een z.g. foutterugkoppel schakeling. Hierop zal 


echter niet nader worden ingegaan. 


In WDF's treden deze effecten niet op indien een over- 
loop nietlineariteit wordt gebruikt met een karakte 


ristiek die ligt in het in fig. 7 aangegeven gebied. 
CONCLUSIES 


Diverse specifieke effecten die kunnen optreden in 

rekursieve digitale filters zijn besproken. We hebben 

gezien dat bij afwezigheid van een ingangssignaal 

oscillaties kunnen optreden die we hebben aangeduid 

met de termen limieteycli en overloop oscillaties. 

Bij excitatie met harmonische signalen trad een sprong- 

verschijnsel op. Bij de realisatie van een digitaal 

filter moeten deze niet gewenste effecten zo mogelijk 

worden vermeden, en diverse oplossingen daarvoor zijn 

aangegeven. Ook al hebben we de genoemde effecten 

echter volledig onder kontrole gekregen dan nog zullen 

de noodzakelijke niet-linesire bewerkingen hun invloed 

blijven doen gelden. Voor kwantisatie resulteert dit 

in kwantisatiefouten, die onder bepaalde omstandigheden 

als (kwantisatie-) ruis kunnen worden beschouwd. Een 

belangrijk aspect bij het ontwerp van een digitaal 

filter is het bepalen van de signaal-kwantisatieruis 

verhouding. Op dit aspect is in het huidige verhaal 

echter niet ingegaan. Voordracht gehouden op 17 februari 1977 op het Natuur- 
kundig Laboratorium Philips, Eindhoven tijdens een 

LITERATUUR gemeenschappelijke vergadering van het NERG (no. 262), 
de Sectie voor Telecommunicatietechniek KIvI en de 

1 H.D. Helms, J.F. Kaiser, L.R. Rabiner (editors), Benelux-Section IEEE. 


"Literature in Digital Signal Processing", IEEE 
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EEN TWEEDE-ORDE DIGITAAL FILTER WAARIN DOOR MIDDEL VAN 
GESTUURDE KWANTISERING "LIMIT CYCLES' VOORKOMEN WORDEN. 


G. VERKROOST 
TECHNISCHE HOGESCHOOL EINDHOVEN 


In this paper a new general second-order digital filter section is proposed. One of the two quantizers 


used in the filter applies controlled quantization, which has the effect that limit cycles are completely 


eliminated under constant-input conditions. 


1. Inleiding. 


In digitale filters kunnen t.g.v. de noodzakelijke 
Kwantisering "limit cycles" optreden (Claasen 1977). 
Dit wil zeggen dat de signalen in het filter na ver- 
loop van tijd niet constant worden indien het. in=- 
gangssignaal constant is. Nu zijn "wave digital fil- 
ters" (Fettweis, 1975) en hieruit afgeleide structuren 
(Meerkötter, 1975) vrij van limit eycles voor een 
ingangssignaal dat gelijk aan nul is, indien "magni- 
tude truncation! wordt toegepast. Anderzijds voorkomt 
“gestuurde kwantisering", in een elementair tweede 
orde filter, limit cycles met een lengte langer dan 
twee voor een willekeurig constant ingangssignaal 
(Butterweck, 1975}, Indien nu gestuurde kwantisering 
wordt toegepast op een bepaalde klasse van wave digi- 
talfilters, worden voor ieder willekeurig constant 
ingangssignaal, limit cycles volkomen uitgesloten 
(Verkroost, 1976), Het is ook mogelijk, uit de voor= 


noemde structuren algemene tweede orde secties af te 


leiden die hetzelfde voordeel hebben (Verkroost, 1977). 


In de onderhavige bijdrage wordt zo'n sectie behan= 


deld, 


In par. 2 wordt het lineaire filter (d.w.z. 
zonder kwantisatoren) geanalyseerd. In par. 3 wordt 
het filter inclusief de kwantisatoren beschouwd. Eerst 
wordt nagegaan hoe, voor het ingangssignaal nul, limit 
cycles kunnen worden uitgesloten. Vervolgens wordt 
aangetoond dat, voor een willekeurig constant ingangs- 
signaal, hetzelfde kan worden bereikt indien gestuurde 
kwantisering wordt toegepast. Om dit artikel beter 
leesbaar te maken, is bij de gegeven afleidingen geen 
gebruik gemaakt van de eigenschappen van de verwante 


“wave digital! structuren. 


2. Analyse van het filter zonder kwantisatoren. 


In Fig. | is het hier beschouwde filter getekend. 

Het ingangssignaal x{(k) veroorzaakt aan de uitgang het 
signaal y{k). Verondersteld wordt dat het ingangssig- 
naal reeds gekwantiseerd is; het kan dus beschreven 


worden met gehele getallen. Het recursieve gedeelte 


van het filter bevat de vermenigvuldigers "Yy Eng: 


Het uitgangssignaal y(k) wordt gevormd door de drie 


Fig. |l Het tweede-orde filter 
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signalen x(k), a, (k) ern -a,(k), na vermenigvuldiging 


pen 
de kwantisatoren aangeduid met 0, en 0, worden de sig- 


met respectievelijk &, , é, en EL op te tellen. Door 


nalen bk) en b, (k) volgens een hierna beschreven 
0 

voorschrift gekwantiseerd, hetgeen de signalen b, Ck) 

resp. b,(k) oplevert. De vertragingselementen zijn 


aangeduid met T. 


We onderzoeken nu eerst het ideale filter 


(zonder kwantisatoren), d.w.z. we beschouwen 0, en 0, 


als doorverbindingen. Uit fig. |. volgt dan 


bk) = a, (ktI) = a,(k) — y‚{ a,(k)ta,(k)tx(k)} (la) 
bk) Er Ta, (k+1) Ee a,(k) gl Yof af | (ke) ta, (lo) tael). (lb) 


Als we de z-getransformeerde van een signaal aanduiden 


met een hoofdletter, volgt uit (1) 
B, (2) = zA, (z) = A, (2) -_ Yí A, CZ)+A, (z)tK(z) | (a) 
B,(z) “=zA,(z) = Aj (z) ok A, CZ)+A, (z)eXlz)}. (2b) 


Met (2) kunnen A, en Á 


5 in X uitgedrukt worden: 


Aj) my, (2) 


RE AFN TEER: 
X(z) z +, 13)2 + Y1*Ya a} 


A, (z) „gert (3a,b) 


eN 
K(z) zZz *lyjergderj trl 


Met (3) volgt dan voor de overdrachtsfunktie H(z): 


Y(z) SjK(z)+8 A, (2)-0,A, (2) 


Hz) = - 


X(z) 


2 
balba jl rf Eolrearel deler 
_Sjz eel (rj-r)-borj bart +6, rytror ddr téjy, 


2 | (4) 
2 +yyradetr,*rg Ì 

We kunnen (4) herschrijven als 
P 

Híz) = cz +dzte (5) 


z“vaztb 
Door een geschikte keuze van Yr Ype 8 8, en 8 3 
kunnen alle gewenste waarden van a t/m e en dus alle 
gewenste overdrachtsfunkties H(z) gerealiseerd worden. 
We beschouwen hier slechts stabiele systemen waarvoor 
moet gelden (zie bijv. Schussler, 1973) 


heat] >o, b-atl zo en b <l. (6) 


Uit deze stabiliteitsvoorwaarden (6) volgt met (4) 


voor de vermenigvuldigers Yy EL Ya: 


Yj >0, Yz >o en Y‚ty, 2. (7) 
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We analyseren nu het recursieve gedeelte van het 
lineaire filter voor het geval dat het ingangssig=- 


naal nul is. 


Met (1) volgt dan dat in de stationaire toestand alle 
signalen in het filter nul zijn. Indien het ingangs- 
signaal echter slechts nul is vanaf een zeker moment 
treden eigentrillingen op. We definieren nu een 


“pseudo-energie' funktie 

Ik) = vb,2(k) + v‚b,2(k) (2) 
en de toename hiervan per tijdstap 

ATl(k) = I(k) — I(k=-1). (8) 


Met (7) zien we dat I(k) niet negatief kan worden: 
verder volgt voor All(k) uit (1) 


AK) = y, (b, 2) - a )} + vlb,l) - a,2(k)}. (10) 


Indien we in (10) met behulp van (1) b,(k) en b,(k) 
uitdrukken in a, (k) en a,(k) (N.B.: x(k) = o) en dit 


uitwerken, volgt 

BIC) = vv Cj*vjD) la, + a, Cm) 
hetgeen met (7) leidt tot 

AT(k) < o. (12) 


De pseudo-energie I(k) kan derhalve niet toenemen in 
de tijd. 


3. Het filter met kwantisatoren. 


We beschouwen nu het filter inclusief de kwantisatoren, 
en wel eerst voor het geval dat het ingangssignaal nul 
is en de kwantisatoren "magnitude truncation' uitvoe- 


ren, hetgeen betekent: 


bol s Ib, | en |b‚GO] < Ib, (|. (13a,b) 
Le) Le 


Dit kwantiseren heeft het effect dat de waarden 
2 
b, (k) en b, 


ter, verkleind worden of tenminste niet vergroot wor- 


Lis, in vergelijking met het ideale fil- 


den; uit (10) blijkt dan dat hetzelfde geldt voor 
A(k). De kwantisatoren werken dus "dissipatief'" op de 
pseudo=energie. Dit betekent dat (12) geldig blijft 
indien de kwantisatoren aanwezig zijn. Het is nog van 
belang om op te merken dat om aan de gelijkheid 

Al(k) = o te voldoen, het tenminste nodig is dat in 
(13a) en (13b) het gelijkteken geldt. Dit wil zeggen 
dat de kwantisatoren dan als doorverbinding werken, 


(k) en b, (k) 
C ü 


hetgeen slechts optreedt indien b, 


gehele getallen zijn. 


We zullen nu bewijzen dat voor x(k) = o geen 
“limit cycles! kunnen optreden. Stel namelijk dat zo'n 
limit eyele wel optreedt, dan bestaat er een getal N 


waarvoor geldt 

b zl di f = ; 
‚{k) oke bk) bk _+N) (l4a,b) 

en dus ook 

Tk) = Tk +N) (15) 


Omdat I(k) niet kan toenemen met toenemende k, kan aan 


(15) slechts voldaan worden indien geldt 
AI(K) = r k=k4 „== 
(k) = o voor k ktl, k +2, k eN, (16) 


hetgeen in ieder geval betekent dat er in dit interval 
niet gekwantiseerd is. Het filter heeft dan dus als 
ideaal filter gewerkt; daar het ideale filter echter 
stabiel is kunnen hierin geen limit cycles optreden. 


De enige mogelijkheid om aan (lá) te voldoen is dan met 
b, (k‚) = b, (k_+N) =O, bk) = b,(k +N)= o, (l7a,b} 


hetgeen wil zeggen dat de stationaire toestand bereikt 
is. De conclusie is dus dat er ondanks de kwantisatoren 
geen limit cycles kunnen optreden in het recursieve 
gedeelte. Het uitgangssignaal y(k) dat een combinatie 
is van de signalen in dit gedeelte is dan ook vrij van 
limit eycles. 

« 
Beschouw nu het filter indien het ingangssignaal x(k) 
de constante waarde X aanneemt. Uit (1) volgt dan voor 


de stationaire toestand 
x(k) = X, a, (k) = b, (k) = =-X en a, (k) = b,(k) = 0.18) 


Om de trillingen in het filter te beschrijven voeren 


we nu verschilvariabelen indie de afwijking tussen het 


werkelijke signaal en dat van de stationaire toestand 


aangeven: 

a, Ck) = a, (k)+x(k) a, (k) = a,(k) (19a,b) 

BK) = b, (l)+x(k) B‚(k) = b,(k) (19e,d) 

Bj Wk) = bj ()tx(k) Bz (k) = bj (k). (19e,E) 
o o o o 


Door substitutie van (19a,b‚e,d) in (la,b}) verkrijgen 
we de vergelijkingen voor deze verschilvariabelen in 


het lineaire geval: 
BU) = a, - yv, la) () + a,()) (202) 


B, (ke) = a, (k) if La, (k) t a, (k)), (20b) 


Indien we nu (Ìl) en (20) vergelijken blijken de ver- 
gelijkingen voor de verschilvariabelen dezelfde te 
zijn als de homogene vergelijkingen voor de oorspron- 
kelijke variabelen (xík) = a). 

Nu is in het voorgaande aangetoond dat bij een in- 
gangssignaal nul, limit eycles worden voorkomen in- 
dien de kwantisering geschiedt volgens het in verge- 
lijking (13) gegeven voorschrift. Hieruit volgt dan 
dat bij een constant ingangssignaal eenzelfde kwanti- 
seringsvoorschrift, nu echter toegepast op de ver- 
schilvariabelen, eveneens limit cycles voorkomt. 

Dit betekent dat de kwantisering dient te geschieden 


volgens 


lk) « [ez |,  (2ta,b) 


le, dl < Ig, O0 
LE, 0 


of, uitgedrukt in de werkelijke signalen: 


b, Ck) + x(k)| slb, (k) + x(k)| 
8 
(22a,b) 


lb,ul s [es (Ke) |. 


Dit houdt in dat in Q, magnitude truncation wordt toec- 
gepast, Het kwantiseringsvoorschrift (22a) voor Q, 
schrijft voor dat Pe wordt gekwantiseerd in de 
richting van (-x(k)), In Q, wordt dus "gestuurde kwan- 
tisering' toegepast, waarbij het ingangssignaal x(k}) 
de kwantisering zodanig stuurt dat indien geldt 

bj, (k) > -x(k) kwantisering "naar beneden! plaats- 
vindt en indien geldt bi) < =x(k) kwantisering 
“naar boven". Indien geldt bj (k) = -x(k), hetgeen 
betekent dat bj, (ke) een zetieel getal is, vindt geen 


kwantisering plaats. 
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Voordracht gehouden op 17 februari 1977 op het 
Natuurkundig Laboratorium Philips te Waalre 
tijdens een gemeenschappelijke vergadering van het 
NERG (n° 262), de Sectie Telecommunicatietechniek 


KIlvIl en de Benelux section IEEE. 


HET LEDENBESTAND VAN HET NERG 


De vitaliteit van een vereniging kan ten dele 
worden afgelezen uit de mate waarin hij appelleert aan 
de belangstelling van de jongeren. Regelmatig melden 
zich nieuwe leden aan en op het eerste gezicht lijkt 
een gunstige diagnose van de levenskracht gewettigd. 
Op grond van een aantal gegevens wordt in het volgende 


dit optimisme enigszins getemperd. 


Aan de hand van de NERG-ledenlijst van augustus 
19/6 en de Naamlijst Ingenieurs 1975 is voor 75% van de 
leden (nl, de 455 die in de Naamlijst zijn vermeld) 
nagegaan in welk jaar zij zijn afgestudeerd, Het 
resultaat ziet u in figuur 1. Het merendeel van de 
academici onder de leden is afgestudeerd in de elektro- 
techniek aan een van de Technische Hogescholen. In 
fig. 2 is het "aanbod" per jaargang van deze categorie 
weergegeven. Als we nagaan welk aandeel per jaargang 
vroeg of laat is toegetreden tot het genootschap dan 


vinden we de grafiek van fig. 3: tussen 1956 en 1968 


schommelt het aandeel rond 10%, daarna treedt een sterke 


daling op. Deze daling kan slechts ten dele worden 
verklaard uit de "incubatietijd" tussen afstuderen en 
toetreden tot het NERG. Uit fig. 4 kan worden afgelezen 
dat de helft van hen die toetreden dit binnen 2h jaar 
doet en driekwart binnen Sk jaar. De sterke daling in 
Fig. 3 zou erop kunnen duiden dat de aanwas van jonge 


leden een stagnatie ondervindt. 


Als de toekomst van het genootschap u ter harte 
gaat kunt u een bijdrage leveren door een wervende 
blik om u heen te werpen en potentiële leden over de 
drempel te helpen. Het is nuttig daarbij te bedenken 
dat het NERG een breed interessegebied heeft dat cok 
vele aanspreekt die een andere studie hebben gevolgd, 


zoals natuurkundigen en geschikte niet-academici. 


G.A. van der Spek 


Secretaris 
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VARIA 


z “Moderne elektronische Meetsystemen'' 


loor de Afdeling voor Electrotech niek van het KIVI 
zal op donderdag I2 en vrijdag 13 januari 1978 te 
Delft een leergang worden gehouden met de titel 
“Moderne elektronische Meetsystemen''. 

Deze leergang wordt georganiseerd in samenwerking 
met de vakgroepen Elektrotechnische Materialen en 
Elektronische Instrumentatie van de Afdeling der 
Elektrotechniek van de TH. Delft, en met medewerking 


van het NERG en de Benelux Sectie IEEE. 


Het doel van de leergang is diegenen, die betrokken 
zijn bij meetinstrumentatie, kennis te laten nemen 
van de recente ontwikkelingen op dit vakgebied. De 
leergang bestaat uit een tiental voordrachten, en 


is onderverdeeld in vier sekties: 


A, Transduktie van informatie 

B, Signaalbewerkingen 

C. Toepassingen van microprocessoren in de 
instrumentatie 


ĳ, Presentatie van informatie 


Het ligt in de bedoeling de komplete teksten van de 
voordrachten te doen uitgeven in de vorm van een boek 
met de tital "Modern Electronic Measuring Systems", 
waarvan de deelnemers van de leergang medio 1978 


een exemplaar zal worden toegezonden. 
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